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1.  INTRODUCTION 

The  objectives  of  this  repon  are  to  present  the  results  of  control  systems  analysis  and  their 
potential  applications  to  control  systems  synthesis  for  high  pcrfomiance  aircraft,  and  identify  a 
proven,  uactablc  and  portable  methodology  to  quantify  the  robust  stability  and  pcrfomiance  of 
highly  integrated  digital  fly-by-wire  control  systems.  These  systems  are  characterized  by  multiple 
distributed  sensing,  processing,  and  actuator  elements  interconnected  by  a  time-division- 
multiplexed  digital  data  bus  layer.  The  aircraft  division  of  Naval  Air  Warfare  Center  at 
Warminster,  PA  is  supponing  this  research  for  developing  an  analytical  methodology  and 
associated  tools  to  assess  the  stability  and  performance  robustness  of  multi-input  multi-output 
(MIMO),  bus-organized  digital  flight  control  systems  The  computation  elements  and  subsystem 
modules  within  the  control  system  may  have  a  multi-rate  structure,  usually  but  not  always,  in  a 
whole  integer  ratio.  Task  synchronization  between  these  elements  may  range  from  no 
synchronization  to  synchronization  at  the  CPU  clock  level. 

This  research  project  consists  of  two  phases  that  are  to  be  performed  sequentially.  The 
results  of  Phase  1  of  the  project  are  presented  in  this  report  which  addresses  the  development  of  a 
robustness  asscssniciii  methodology  suitable  for  analyzing  a  generic  class  of  MIMO,  bus- 
organized  flight  control  systems.  The  future  work  in  Phase  2  is  projected  to  deal  with 
implementation  and  validation  of  the  methodology  developed  in  Phase  1.  This  implementation  in 
Phase  2  should  account  for  bus  data  delays,  multiple  distributed  computation  delays,  and  various 
degrees  of  synchronism  (or  lack  thercoQ  between  the  distributed  processing  elements  and  the  data 
bus.  This  report  is  intended  to  provide  useful  infomiaiion  to  the  engineers  and  scientists  from 
different  organizations  of  DoD,  NASA,  and  leading  aerospace  companies. 

The  currently  available  theories  and  methodologies  for  robustness  assessment  of  multi-input, 
multi-output  (MIMO)  systems  have  been  reveiwed  and  thoroughly  examined  in  this  repon.  Recent 
literature  has  shown  three  major  different  approaches  to  address  the  robust  stability  and 
perfomiance  of  MIMO  control  systems.  The  first  approach  relies  on  the  structured  singular  value 
(p)  and  Hoo  analyses  in  the  framework  of  Linear  Fractional  Transformations  (LFT)  of  MIMO 
systems  that  are  subjected  to  uncenainties  with  known  bounds  [Doyle  (1982)]  and  [Safonov 
(1982)].  The  concepts  and  fundamental  results  and  applications  of  ^.-analysis  and  p-synthesis  arc 
presented  in  the  publications  of  Doyle  and  his  coworkers  (Skogestad  et  al.  (1988),  Stein  and  Doyle 
(1991),  Balas  et  al.  (1991),  Packard  and  Doyle  (1993),  Packard  et  al.  (1993)].  The  second 
approach  [Yeh  et  al.  (1985a),  Molander  and  Willems  (1980),  Tao  et  al.  (1991)]  generalizes  the 
definition  and  calculation  of  gain  and  phase  margins  for  MIMO  systems  analogous  to  those  for 
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single-input,  single-output  (SISO)  systems.  The  third  approach  jFreudenberg  and  Looze  (19S(i), 
Freudcnbcrg  (1990)]  examines  the  relationship  between  the  open  loop  gain  and  phase  with  respect 
to  the  closed  loop  sensitivity  and  complementary  sensitivity  functions  lot  SISO  systems,  and  then 
formulates  heuristic  rules  for  shaping  the  loop  transfer  functions  for  MIMO  systems.  We  have 
emphasized  the  first  approach  of  |i-analysis  for  robustness  assessment  of  MIMO  systems,  sshich 
constitutes  the  main  txxly  of  this  report.  A  simulation  example  has  been  provided  to  illustrate  its 
usage.  However,  technical  discussions  on  the  second  and  third  approaches  are  presented  in  the 
appendices  without  any  simuFtion  examples. 

The  stability  robustness  tests  for  MIMO  systems  stans  with  the  basic  requirement  that  the 
nominal  closed-loop  system  must  be  stable,  which  can  be  tested  via  the  multivariable  Nyquist 
stability  theory.  Furthermore,  stability  robustness  tests  are  also  based  on  the  multivariable  Nyquist 
theory  which  is  briefly  described  in  one  of  the  appendices. 

This  report  is  organized  in  five  sections  and  three  appendices.  Section  2  presents  the  basic 
concept  of  the  structured  singular  value  (p)  and  its  properties  along  with  the  techniques  of  p- 
analysis  and  p-synthesis.  The  results  of  simulation  experiments  for  the  flight  control  system  of  an 
advanced  aircraft  are  discussed  in  Section  3  to  illustrate  the  efficacy  of  the  p-analysis  technique  for 
MIMO  systems.  Section  4  desenbes  the  codes,  developed  in  the  MATLAB  environment,  which 
can  be  used  for  robustness  analysis  and  synthesis  of  MIMO  control  systems.  The  last  section 
summarizes  and  concludes  the  findings  of  this  research  report.  Appendix  A  briefly  describes  the 
principle  of  multivariable  Nyquist  criterion  which  is  essential  for  understanding  the  stability  of 
closed  loop  control  systems.  Appiendix  B  provides  a  general  definition  and  computational  methods 
for  gain  and  phase  margins  of  MIMO  systems  in  a  framework  analogous  to  those  for  SISO 
systems.  Appendix  C  describes  the  relationship  between  open  loop  gain  and  phase  with  respect  to 
the  closed  loop  sensitivity  and  complementary  sensitivity  functions  for  SISO  systems,  and  then 
formulates  heuristic  rules  for  shaping  the  loop  transfer  functions  for  MIMO  systems. 
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2.  STRUCTURED  SINGULAR  VALUE  (n)  ANALYSIS 

This  section  introduces  the  technique  of  robustness  analysis  for  MIMO  control  systems 
based  on  Hoc  and  ^.-analysis  &  synthesis, 

2.1.  The  Structured  Singular  Value 

Figure  2-1  shows  how  the  plant  pevturbations  represented  by  A(s),  .ve  C,  interact  with  the 
finite-dimensional,  linear,  time-invariant  control  system  M(s)  which  is  the  based  on  the 

nominal  plant  model.  The  input  t//  to  the  control  system  M  consists  of:  all  exogenous  signals, 
namely,  the  reference  command(s)  to  be  tracked,  disturbances  and  sensor  noise;  and  the  feedback 
control  input.  The  output  ^  of  the  control  system  M  consists  of;  all  plant  variables  needed  for 
specifying  the  stability  and  performance  criteria;  and  the  sensor  data  feeding  the  controller.  In  the 
definition  of  the  structured  singular  value  jJ.^{M(s^))  of  the  transfer  matrix  M(s)  at  a  given  Sq  , 

the  underlying  uncertainty  A(s)  belongs  to  a  set  of  matrices,  A{s},  which  is  prescribed  to  have  a 
block  diagonal  structure  with  the  following  three  characteristics: 

•  Type  of  each  block; 

•  Total  number  of  blocks; 

•  Dimension  of  each  block. 

In  general,  there  are  two  types  or  blocks;  Repeated  scalar  blocks  and  full  blocks.  Let  two 
nonnegative  integers,  S  and  F,  represent  the  number  of  repeated  scalar  blocks  and  the  number  of 
full  blocks,  respectively.  We  introduce  positive  integers,  to  represent  the 

dimensions  of  these  blocks  such  that 

•  The  1^  repeated  scalar  block  is  3,!^^  where  is  the  x  identity  matrix  and  5^  e  C 

•  The  full  block  belongs  to  . 


Figure  2-1.  The  Closed-Loop  Control  system 


Definition  I:  The  block  diagonal  structure,  A(5)  ,  is  defined  as: 

4(.S)  =  ^dia}>^5\(s)lr^  ,...,Ssts)lr^,A]{s) . A/r(5)|:6,  eC,Aje  |  c 


(1) 
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The  following  constraint  must  be  satisfied  for  consistency  of  the  dimensions: 


5  F 

X''' 

,=1  ;=1 

The  following  notation  is  introduced  in  'erms  of  norm-bounded  subsets  of  A  : 


BA=  {A  e  A\a(A)  <  1} 


(3) 


Definition  2:  For  anyM  e  C"^",  its  structured  singular  value  is  defined  as; 


_ \ _ 

inf|a(d);A  e  A,dci{I  -  MA)-  0; 


(4) 


[0  Vd  e4,det(/-Md)^0 

Remark  I ;  The  physical  significance  of  fj.^  (M)  can  be  interpreted  following  the  feedback  loop  in 
Figure  21  where  the  loop  equations  are: 

<^  =  A/i/r;  and  y/  =  A^  (5) 


As  long  as  I-MA  is  nonsingular,  the  only  solution  for  the  loop  equations  is  y/  =  ^  =  0. 
However,  if  /  -  MA  is  singular,  then  there  are  infinitely  many  solutions,  and  the  norms  II  y/W  and 

li  ^  II  can  be  arbitrarily  large.  Thus  the  feedback  sy.stem  becomes  unstable.  Hence,  (M))  ^  is 
a  measure  of  the  smallest  A  belonging  to  a  given  class  of  uncertainty  stiucture  that  causes 
instability  of  the  above  MIMO  feedback  loop.  In  other  words,  (ki^(M))  is  the  size  of  the 
smallest  "destabilizing  "  A  q  A  which  satisfies  the  condition  of  det(l  -  MA)  =  0.  It  follows  from 
the  MIMO  Nyquist  stability  test  (see  Appendix  A)  that  the  zeros  of  dei{l  -  MA)  are  the  closed- 
loop  poles  of  the  feedback  system.  Therefore,  if  M  =  Af(5())  is  a  transfer  matrix  evaluated  at  a 

point  sq  in  the  complex  plane,  (fd^  (M ))  *  is  the  size  of  the  smallest  allowable  A  which  moves  a 
closed-loop  pole  to  that  location  sq. 


The  structured  singular  value  id^{M)  of  a  transfer  matrix  M  can  be  related  to  the  familiar 
linear  algebraic  quantities  if  A{s)  is  bounded  for  every  seC  as  seen  below  (Doyle  et  al.  (1982)]: 

Result  J:  If  A  -{5l:5eC}(S  =  ],F  =  0,rj=n).  then  (M)  =  p(M)  where  p(M)  is  the 
spectral  radius  of  M. 
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Result  2:  If  A  (S  =  0,F  =  l,mj  =n),  then  o{M)  where  o{M)  is  the  largest 

singular  value  of  M. 

Result  3:  For  a  general  A,  the  following  inequality  holds. 

{5l„:5eC]^AcC''^'‘  (6) 

From  the  definition  of  fi ,  and  Results  1  and  2  above,  u  follows  that: 

p(M)<^^(M)<cf(M)  (7) 

vSince  the  difference  between  p  and  c  can  be  very  large,  the  bounds  in  Result  3  may  not  be 
sufficiently  light  for  evaluation  of  p^(Af).  The  following  result  [Doyle  (1982)]  snows  how  these 

bounds  can  be  refined  by  transformations  on  M  that  do  not  affect  p^{M)  but  do  affect  p  and  o . 
Result 4:  The  following  propeny  of  holds: 


=  ^a{QM)  -  p^{M)  -  p^[dMD~^) 

(8) 

with  Q  €Q 

QAeA  o(QA)  =  ^(AQ)  =  a(A)  DA  =  AD 

(9) 

for  any  A  A  . 

,  D  &  D  and  QeQ  where  the  sets  D  and  Q  are  defined  as: 

D  =  |d(ag[c>, 

■  DsA'r,, . ]  A  6  =  £)•  >  O.dj  e  R.dj  >  o) 

q  =  {q 

(10) 

FurJiermore  the  bounds  of  (M)  in  (7)  can  be  tightened  to 
max  max  inf 

C  .  4  ^  ^  -  D  .  )  <"> 

It  is  shown  that  the  lower  bound  in  eq.  (1 1)  is  an  equality  [Doyle  (1982)].  Apparently,  no 
analytical  techniques  exist  to  find  the  global  maximum  of  the  spectral  radius  p(QM)  although  it 
may  have  local  maxima.  Thus  local  seaich  cannot  be  guaranteed  to  obtain  p  from  its  lower  bound. 
The  p  software  package  [Balas  et  al.  (1991))  uses  a  slightly  different  algorithm  to  compute  the 
lower  bound.  While  th„re  arc  open  questions  about  its  convergence,  this  algorithm  usually  works 
well  and  has  proven  to  be  an  effective  tool  for  approximate  computation  of  p. 

The  upper  bound  can  be  reformulated  as  a  convex  optimization  problem  so  that  the  global 
minimum  can  be  found  in  principle.  Unfortunately,  the  least  upper  bound  is  not  always  equal  to 


5 


NAWCADWAR-93010-60 


|i.  For  block  structures  4  satisfying  2S  +  F  <  3,  the  least  upper  bound  has  been  shown  to  be 
always  equal  to  and  for  block  structures  with  2S  +  F  >  3,  |i  may  not  attain  the  least 

upper  bound.  These  results  are  summarized  in  the  Table  2-1  to  show  for  which  cases  the  least 
upp>er  bound  is  guaranteed  to  be  equal  to  p. 


Table  2-1.  Conditions  for  p  reaching  its  upper  bound 


0 

1 

2 

3 

m 

0 

— 

yes 

yes 

yes 

no 

1 

yes 

yes 

no 

no 

no 

2 

no 

no 
— — 

no 

no 

no 

!t  is  essential  to  estimate  the  upper  and  lower  bounds  for  reliable  use  of  the  p  theory.  The 
most  important  usage  of  the  upper  bound  is  in  a  computational  scheme  when  it  is  combined  with 
the  lower  bound.  The  other  important  feature  of  the  upper  bound  is  that  it  can  be  combined  with 
//oo  control  synthesi.s  methods  to  yield  an  iterative  p-synthe.sis  method. 

2.2.  Linear  Fractional  Transformation 

'i  he  dynamic  system  model  in  Figure  2-2  includes  both  parametric  and  non-parametric 
uncertainties  where  the  complex  matrix  M  of  the  nominal  system  is  partitioned  as: 


Mjj 

M21 


Mi2 

^22 


(12) 


and  the  uncenainty  block  stnicture  A  is  compatible  in  size  with  M22-  For  4  loop 

equations  are: 

z  =  Mi\d -h  Mi2  p.'  w  =  Mjxd  +  MxiP:  and  p  =  Aw  (13) 


Figure  2-2.  Linear  Fractional  Transformation  Structure 
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Definition  3:  The  system  in  Figure  2-2  is  well-posed  if  there  exist  unique  vectors  z,  w,  and  p 
such  that  the  equation  set  (13)  is  satisfied  for  any  exogeneous  input  vector  d.  In  other  words,  the 
system  is  well-posed  if  and  only  if  /  -  M22^  is  invertible. 


Remark:  2:  If  the  system  in  Figure  2-2  is  not  well-posed,  then  there  are  infinitely  many  solutions 
to  the  loop  equation  set  (13). 


Definition  4\:  If  the  system  in  Figure  2-2  is  well-posed,  then  the  linear  fractional  transformation 
(LFT),  Ff(M,A),oi\M  induced  by  ded  is  defined  as: 

Ff  (M ,A )  —  ( ^  ~  ^21^ )  ^2\ 

and  (14) 

z  =  F((M.A )d  is  satisfied. 


Remarks:  The  subscript  /  on  Ff  pertains  to  the  "lower "  loop  of  M  being  closed  by  d  e  d  as 
shown  in  Figure  2-2.  An  analogous  formula  describes  F^^(M,A),  which  is  the  resulting  matrix 
obtained  by  closing  the  "upper"  loop  of  M  via  an  uncertainty  matrix  having  z  as  the  input  and  d  as 
the  output. 


Remark  4 :  In  the  LFT  formulation  of  eq.  (14),  the  sub-matrix  Mjj  represents  the  nominal  plant 
transfer  matrix;  d  e  Bd  ,  as  defined  in  eq.  (3),  can  be  viewed  as  a  norm-bounded  perturbation 
from  an  allowable  set;  and  the  sub-matrices  M]2,  M2]  and  M22  provide  the  information  of  how 
the  parametric  and  non-parametric  uncenainties  influence  Mjj. 


Given  two  defined  block  structures,  and  A 2,  which  are  compatible  in  size  with  Mj j  and 
M22  respectively,  a  block  structure  A  is  generated  as: 


t  d;^,d2  e  A2 


(15) 


Now  the  structured  singular  value  p  can  be  computed  relative  to  the  above  three  block 
structures,  d,  dj  and  do.  In  general,  the  following  two  problems  are  of  interest: 

•  Uetermination  of  whether  the  LFT  is  well  posed  for  all  A2  e  A  2  with  <7(d  7 )  ^  i  :  and 

•  Evaluation  of  a  scalar  measure  of  die  transfer  matrix  Ff(M,A2)  relative  to  the  nonn-boundcd 
set  of  penurbations  4.2. 

To  solve  the  above  two  problems  by  calculating  the  structured  singular  values  (p),  we  present  two 
theorems  ([Packard  et  al.  (1993)]  and  Balas  et  al.  (1991)])  which  form  the  foundation  of 
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^-analysis  and  synthesis.  In  this  respect,  the  following  notations  are  used;  H](*)  is  with  respect 
to  4.1.  to  ^2'  ^tid  pi;  if*)  is  with  respect  to  4  ss  defined  in  eq.  (15). 

Theorem  1  (Theorem  3.4,  Balas  et  al.  (1991),  p.  2-13):  TheLFT  Ff(M,A2)  is  well  posed  for  all 
^2  ^  ^A2  if  ottly  if  P2  (^22  ) 

Theorem  2  (Theorem  3.5,  Balas  et  al.  (1991),  p.  2-13);  if  and  only  the  following 

two  conditions  hold: 


•  ^^2(!^22)<  J 

Max  Hj (F ^(M ,A2))  <  1 

A2  e  BA2 


(17) 


2.3.  Review  of  p.  in  the  Frequency  Domain 

The  most  well  known  use  of  p  as  a  robustness  analysis  tool  is  in  the  frequency  domain.  Let 
M(s)  be  the  nominal  transfer  matrix  of  a  stable,  multi-input,  multi-output  transfer  function  of  a 
linear  time-invariant  control  system.  Given  the  block  structure  of  all  penurbations  which  are 
themselves  dynamic  systems  having  the  structure  of  A,  the  problem  is  to  determine  the  effects  of 

these  perturbations  on  the  control  system.  To  accomplish  this  task,  first  let  E  denote  the  set  of  all 
real-rational,  proper,  stable,  transfer  matrices.  Associated  with  any  block  structure  A,  let  0{A) 

denote  the  set  of  all  block  diagonal,  stable  real-rational  transfer  matrices,  with  block  structure  of 
2I. 


0(4)  =  {4(^)€  5.-  jeC+} 


(18) 


.  _  ^ 

Uncertainties 

A(s) 

r 

Nominal  System 
M(s) 

^ _ 

Figure  2-  3.  Linear  Fractional  Transformation  in  the  Frequency  Domain 

3  (Theorem  3  6,  Balas  et  al.  (1991),  p.  2-14):  For  p>0,  the  feedback  loop  shown  in 
Figure  2-3  is  well-po.sed  and  internally  stable  for  all  A  e  0(4)  with  ||4||^  <  ft  ^  if  and  only  if 
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^  ^A^(jco))<f}  (19) 

Remark  .5::  Theorem  3  ensures  th?.t  the  peak  value  on  the  ji-plot  of  the  frequency  response 
corresponds  to  the  penurbation  whose  size  determines  the  robust  stability  of  the  loop  system. 


Often  performance,  in  addition  to  stability  of  a  closed-loop  system,  must  be  robust  to 
perturbations.  Typically  there  are  exogenous  disturbances  acting  on  the  system  which  result  in 
tracking  or  regulating  errors.  Under  the  penurbation,  the  effects  of  these  disturbances  on  the  error 
signals  can  significantly  increase  and  consequently  the  closed-loop  system  pierformance  may  not  be 
acceptable.  Therefore,  a  procedure  for  robust  performance  test  is  necessary’.  Such  a  test  should 

indicate  the  worst  case  performance  degradation  associated  with  a  given  level  of  perturbations.  If 
the  transfer  matrix  M  is  stable,  real-rational  and  proper  with  rip  inputs,  and  n^.  +  outputs, 

then  M  can  be  partitioned  such  that  Mjj  has  ttp  inputs  and  outputs,  M22  has  inputs  and 

outputs,  and  the  remaining  two  submatrices  of  M  have  compatible  dimensions.  Let  the  block 

structure  of  the  uncertainties  corresponding  to  the  system  submatrix  Mjj  be 

augmented  as  follows: 


~sys 


-1 

\Aj 

01 

V 

i - 

0 

1 _ 

“  ^•42  cC"'' 


The  above  block  structure  is  set  to  address  the  robust  perfomiance  of  the  loop  shown  in  Figure 
2-4,  and  the  perturbed  (closed  loop)  transfer  matrix  is  denoted  by  F^(M,d)  with  A(s)e 


Theorem  4:  (Theorem  3.7,  Balas  et  al.  (1991),  p.  2-161:  Let  >  0.  For  all  zl(5)  e  with 

II^IL  ^he  loop  in  Figure  2-4  is  well-posed,  internally  stable,  and  ||F„(Af  ,d)||^  <P  if  and 

H 

only  if 


(21) 


In  general,  a  procedure  for  p-analysis  and  p-synthesis  of  a  control  system  to  assure  robust 
performance  design  procedure  involves  several  steps: 

Step  1 :  Specification  of  the  closed  loop  feedback  structure. 

Step  2:  Specification  of  the  modeling  uncertainly  and  performance  objectives  in  terms  of 
freque^'cy-dependent  weighting  matrices. 

Step  3:  Construction  of  open-loop  interconnection  for  the  control  synthesis  routines. 

Step  4:  Design  of  the  connoller  for  the  open-loop  interconnection  in  Step  3. 
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Step  5:  Analysis  of  robustness  properties  of  the  resulting  closed-loop  systems  using  the  structured 
singular  value,  known  as  p-analysis. 

Step  6:  Identification  of  frequency-dependent  matrices,  obtained  in  the  p-analysis  step,  to  scale 
the  open  loop  interconnection. 

Step  7 :  Rede  sign  of  the  scaled  Hoo  controller. 

Step  8:  Iteration,  if  necessary,  on  Steps  5,  6  and  7,  known  as  p-synthesis. 

d  Disturbance 
p  Perturbation 
u  Control 

w  Weighted  control 
z  Performance  variable 


Figure  2-4,  Transfer  Matrix  with  Penurbations 

2.4,  The  p— Analysis: 

The  objective  is  to  analyze  a  stabilizing  controller  K  for  all  stable  perturbations 
with  <  7.  The  controller  K  must  satisfy  the  following  conditions: 

•  The  perturbed  closed-loop  system  remains  stable 

•  The  perturbed  weighted  sensitivity  transfer  matnx  relating  z  \o  d 

S(AG)^Wpll  ^G(/  ■hAcW^^l)Kr^  (22) 

has  <  7  for  all  such  perturbations. 
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Uncertainty  Models:  Given  the  nominal  plant  transfer  matrix  G ,  it  is  necessar>'  to  specify  an 
uncertainty  weighting  matrix,  Wdel .  and  a  normalized  variable  transfer  matrix  Both 
and  Aq  have  to  be  stable  transfer  matrices,  and  there  must  not  be  any  pole-zero  cancellation 
between  ,  Aq  and  G  in  the  closed  right-half  plane.  All  of  the  uncenainties  in  modeling  the 
plant  should  be  captured  in  the  normalized,  variable  transfer  matrix  Aq,  which  is  used  to 
parametrize  the  postulated  differences  between  the  nominal  model,  G ,  and  the  actual  behavior  of 
the  real  plant  denoted  by  G  =  G(l  +  AqWj^i )  Therefore,  the  two  transfer  matrices  ^del  Aq 
parametrize  an  entire  set  of  plants  which  must  be  appropriately  controlled  by  the  robust  controller 
K.  The  actual  plant  transfer  matrix  which  includes  the  modeling  errors  and  disturbances  belongs  to 
the  following  set; 

X  =  {G(I  +  AgW^i}:Ag  stable,  \\AgL  (23) 


The  performance  of  the  closed  loop  nominal  system  is  evaluated  using  the  sensitivity  matrix, 
(I  +  GK)~^  between  the  disturbance  d  and  the  performance  variable  z  in  Figure  2-4.  The  objective 
is  to  chr  enze  good  performance  in  terms  of  a  weighted  -norm  on  this  transfer  manix. 
Given  a  stable,  real-rational  transfer  matrix  VVp,  the  nominal  performance  is  achieved  if 


\V p(I  +GK)  '  <1.  Similar  to  the  uncenainty  model  where  the  frequency-dependent 


weighting  functionW^e/  normalizes  the  specifications  Aq  for  plant  modeling  unci  taintics,  the 
frequency-dependent  weighting  function  Wp  is  used  to  normalize  the  performance  'Uch  that  the 
norm  of  the  weighted  sensitivity  matrix  is  less  than  1 . 


To  analyze  the  controller  of  a  given  plant  which  is  subjected  to  uncenainties,  four  different 
enteria  need  to  be  tested  as  explained  below: 


•  Nominal  Staoilitv:  The  control  system  with  the  nominal  plant  in  the  loop  is  internally  stable. 
That  is,  there  are  no  pole- zero  cancellations  in  the  closed  right-half  r-plane,  and  the  poles  of 
the  closed-loop  system  transfer  matrix  must  lie  in  the  open  left-half  5-plane. 


•  Nominal  Performance:  The  performance  objective  is  satisfied  for  the  nominal  plant  model  if: 


Wp(I+GK)-^ 


<  I 


(24) 


•  Robust  Stability:  The  closed  loop  system  between  the  penurbation  p  and  the  uncertainty 
weighted  control  vv  -Wdel  ^  internally  stable  for  all  possible  plant  models  Gel,  i.e., 
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\^delf^G(I  +  KG)-^\^<1  (25) 

Robust  Performance:  The  closed-loop  system  is  internally  stable  for  all  G  e  £  and,  in 


addition,  the  performance  objective  is 

W  {I  ^GKr^\  <1  VGel  (26) 

I  oo 

Remark  3:  Referring  to  Figure  2-4,  nominal  performance  can  be  equivalently  expressed  as: 

W^llL  <  1  (27) 

and  robust  stability  as: 

\\^^(^1I)\L<^■  (28) 

The  following  are  the  conditions  of  robust  stability  and  performance; 

11miCMii)|L<1  and  (29) 

which,  by  Theorem  2,  are  equivalent  to 

(30) 


2.5.  The  4-Synthesis 

It  follows  from  the  4  analysis  that  the  inequality  4j  <  a(DA/D~*)  holds  for  any 
diagonal  matrix  D  with  a  block  structure  matching  that  of  A  =  diaj?(Ai,...,A„).  In  the  4-synthesis 
procedure  [Packard  and  Doyle  (1993)],  the  strategy  is  to  pick  stable  diagonal  transfer  matrices  D(s) 
with  stable  inverses,  then  to  use  //«,  synthesis  to  compute  the  controller  K  which  minimizes 
||Df^(G,^;D"^||  .  TheD  -matrices  may  be  either  constants,  with  no  state  dynamics,  or 

frequency  dependent.  The  latter  case  would  obviously  increase  the  number  of  states  in  the 
synthesized  controller  K,  and  thereby  add  their  states  to  those  of  the  closed  loop  control  system. 

The  4-synthesis  procedure  can  be  executed,  after  the  initial  con'  -vl’er  design,  as  a  two-phase 
iterative  process  called  the  D-K  iteration  [Packard  et  al.  (1993)]  which  is  summarized  below: 
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Phase  0;  Design  a  conooller  K(s)  based  on  a  given  nominal  plant  transfer  matrix  (which  may 
include  the  performance  and  stability  weighting  functions)  such  that  the  closed  loop 
system  is  internally  stable  and  satisfies  the  criteria  of  performance  specifications. 

Phase  1:  Execute  the  p-analysis  procedure  using  the  closed  loop  system.  If  the  robu.stness 
requirement  is  not  satisfactory,  then  pick  a  diagonal  D(s)  so  that  D(s)  andD~'(i')  are 
stable,  and  the  upper  bound  <j( DFi(G,K }D~^ )  is  approximately  minimized  ever  the 
frequency  range  of  interest. 

Phase  2:  Execute  synthesis  to  identify  a  controller  K(s)  which  minimizes  the  norm 

A' of  the  scaled  system.  Go  back  to  Phase  1  and  repeat  the  iterative 

procedure  until  the  desired  condition  of  <  1  is  achieved  over  the 

frequency  range  of  interest  or  no  further  reduction  in  p  is  possible. 
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3.  A  SIMULATION  EXAMPLE 

The  (continuous-time)  flight  dynamic  model  (Safonov  et  al.  (1981)]  of  an  advanced  aircraft, 
linearized  at  tlie  operating  condition  of  7.62  kilometers  and  0.9  Mach,  was  selected  to  illustrate  the 
ja  analysis  and  ^.-synthesis  procedures.  The  plant  matrices  are  given  below. 

■-O.0226  -36.6170  -18.8970  -32.0900  3.2509  -0.7626" 

0.0001  -1.8997  0.9831  -0.0007  -0.1708  -0.0050 

0.0123  11.7200  -2.6316  0.0009  -31.6040  22.3960 

A  = 

0  0  1.000  0  0  0 

0  0  0  0  -30.0000  0 

0  0  0  0  0  -30.0000 

"0  1  0  0  0  0" 

0  0  0  10  0 

The  six  plant  state  variables  are  forward  speed,  angle  of  attack,  pitch  rate,  attitude  angle,  eleven 
actuator  position,  and  canard  actuator  position;  the  two  control  inputs  are  eleven  and  canard 
signals;  and  the  two  output  variables  are  angle  of  attack  and  atthude  angle.  The  open  loop  plant 
poles  are  given  in  the  Table  3-1..  Therefore,  the  nominal  plant  is  unstable  in  the  open  loop  with  a 
pair  of  complex  eigenvalues  in  the  right-half  5-plane. 

Table  3-1.  Plant  open-loop  poles 
Open-loop  poles  of  the  plant 


-5.6754 

6.9{)02e-l-t-2.48e-li 

6.9002e-l-2.48e-ii 

-2.5794e-l 

-3e-t-l 

-3e-Hl 


0  0  0  0  0  30 
0  0  0  0  30  0 


A  diagram  of  the  closed-loop  system,  which  includes  the  feedback  structure  of  the  plant  and 
the  controller,  and  elements  associated  with  the  uncertainty  models  and  performance  objectives,  is 
shown  in  Figure  3- 1 . 

The  dashed  box  represents  the  actual  plant,  with  associated  transfer  function  G.  Inside  the 
box  is  the  nominal  model  of  the  airplane  dynamics,  G ,  and  two  elements,  and  Aq  ,  which 
represent  the  uncertainties  in  the  plant  model.  The  frequency  weighting  function  is  assumed 
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to  be  known  and  reflects  the  amount  of  uncertainty  in  the  model.  The  transfer  matrix  is 
assumed  to  be  stable  and  unknown,  except  for  the  norm  condition,  ||/^g'IL  performance 

objective  is  that  the  transfer  matrix  from  d  to  z  be  small,  in  the  |1*11„  sense,  for  all  possible 
uncertainty  transfer  matrices  Aq  .  The  weighting  function  reflects  the  relative  irnponance  of 

various  frequency  ranges  for  which  the  performance  is  desired. 


Figure  3-1.  Closed-Loop  System  of  the  Aircraft  Control  with  Uncertainties 

The  control  objective  is  to  identifv  a  stabilizing  controller  K  such  that  for  all  stable 
perturbations  Aq{s),  with  ||d(;|^  <  1,  the  penurbed  closed-loop  system  remains  stable,  and  the 

perturbed  weighted  sensitivity  transfer  matrix, 

S(Ac)  =  Wp(l  -h  6(1  +  AGW^i)Kr^ 

has  ||5(A(;)j|^  <  1.  It  follows  from.  Section  2  that  these  objective  functions  are  compatible  to  the 
structured  singular  value  framework. 

To  formulate  the  problem  into  the  framework  of  (J-analysis,  we  augment  the  perturbed  plant 
as  the  oi>en-loop  system  structure  of  p_block  shown  in  Figure  3-2. 

At  the  initial  stage,  we  used  three  conventional  algorithms  for  multivaiiable  control  systems 
synthesis,  namely,  LQG,  H2,  and  Hoo.  Our  intention  was  to  compare  the  robustness  of  these  three 
controllers  under  the  same  structure  and  bound  of  uncertainty  blocks  and  disturbance  blocks,  and 
then  formulate  a  strategy  for  quantitative  evaluation  of  robust  stability  and  performance.  The  results 
are  summarized  below. 
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weighted 


perturbation  p 
disturbance  d 
control  u 


weighted  control  w 
performance  variable  z 
error  signal  e 


Figure  3-2.  Augmented  Open-Loop  System 

3.1.  The  LQG  Controller 

In  the  LQG  synthesis,  it  is  well  known  that  the  nominal  closed-loop  system  is  stable  when 
controlled  by  the  LQG  compensato.*.  The  closed-loop  poles  of  the  nominal  system  are  given  in 
Table  3-2  for  the  state  weighting  matrix  Qc=C3^C;  control  weighting  matrix  Rc=l2;  plant  noise 
covariance  matrix  QpjBB^^;  and  measurement  noise  covariance  matrix  Rf=10‘^  l2.  As  the  nominal 
plant  model  was  augmented  with  the  modeling  uncertainty  blocks,  the  resulting  closed  loop  system 
under  LQG  control  became  unstable.  The  closed  loop  poles  are  shown  in  Table  3-3.  This 
instability  results  from  the  uncertainties  imposed  on  the  nominal  plant.  In  that  case,  the  loop 
transfer  matrix  from  perturbation  p  to  weighted  control  w  in  Figure  3-2  (after  closing  the  loop  with 
the  LQG  compensator)  is: 

w  =  Kis}(l  -G{s)K{s))~^G(s)p 

The  LQG  system  does  not  guarantee  that  the  closed-loop  system  will  remain  stable  in  the  presence 
of  uncertaintes  and  disturbances  [Doyle  (1978)1  although  LQR  guarantees  a  gain  margin  of  (•j,«>) 

and  phase  margin  of  (-60^,-^60'’).  In  this  particular  case,  the  output  becomes  easily  unbounded 
because  G{s)  is  unstable  in  the  open  loop. 
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Table  3-2.  Closed-loop  poles  by  Table  3-3.  Closed  loop  poles  by 
LQG  Design  with  by  LQG  Design  witii 

Exact  Parameters  Uncenain  Parameters 


Closed-loop  poles  by  LQG  design 

Closed-loop  poles 
of  augmented  structure 

MA-BF) 

X{A-LC) 

-5.8348+4.4879i 

-5.8348-4.4879i 

-1.3878 

-3.0056e+l 

-3.00c-t-l 

-20.199-H28.714i 

-20.199-28.714i 

-4.2003e-t-l 

-2.1404e-2 

-2.6961 

-2.99l0e-H 

-2 1.444-1-37.6991 
-2 1.444-37.6991 
-5.2851e+l 

7.3409 

-6.3281 

-2.7026 

-2.9909e-t-l 

-3.00e-t-l 

-8.271e-l 

-2.9987e-t-l 

7.0703e-3 

-2.1401e-2 

3.2.  H2  and  //»  Controller  Analysis  Using  }x 

Two  controllers  were  synthesized  via  and  //_  norm  optimization  to  stabilize  the 

augmented  open-loop  structure  in  Figure  3-2.  Based  on  the  plant  model  and  the  desired 
performance  of  the  closed  loop  control  system,  the  respective  weighting  matrices,  and 

for  uncertainty  and  performance  were  chosen  to  be: 


J  -4-  3 

0 

s  +  3 

■  ; 

0 

1 

II 

s  + 300 

0 

and  V*/p  = 

s  +  0.03 

0 

s  +  300. 

s  0.03 . 

These  two  controllers  were  compared  using  the  p-analysis  criteria  described  in  Section  2.  The 
selected  set,  Sp.of  frequency  points  at  which  the  control  systems  were  analyzed  and  synthesized 
consists  of  f  -  60  points  logarithmically  equally  spaced  between  the  decades  of  10  to  10"*.  The 
results  for  three  different  cases  are  shown  in  Tables  3-4,  3-5,  and  3-6,  respectively.  The 
uncertainty  and  the  external  disturbam  es  in  Case  1  are  represented  by  two  full  blocks,  i.e., 
4)  and  4.2  This  imj-lics  that  the  uncertain  parameters  in  the  plant  could  be 

perturbed  arbitrarily  in  ariy  form  within  a  given  bound.  The  parameter  perturbation  in  the  designed 
compensator  is  analogous  to  the  concept  of  gain  margin  and  phase  margin  in  the  SISO  systems. 
As  seen  in  Table  3-4,  both  H.  and  controller:'  do  not  show'  good  nominal  and  robust 
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perfonnance  (||3/22lL  ‘ire  larger  than  1)  aithough  the  controller  is  superior 

J<i<i 

to  the  //j  controller  from  these  perspectives.  However,  the  stability  robustness  is  good  for  both 
controllers  because  max  ^(Mii((Oi))  is  smaller  than  1  for  each. 

i<i<i 


Table  3-4.  Comparison  of  the  controllers  (Case  I) 


Controller 

Algorithm 

robust  statibiliity 

robust  performance 

IW22L 

l<i<t 

l^/L 

max  n(M((i)i)) 
I<i<l 

l|4L 

■nn^ym 

1.5154 

3.0918 

1.8184 

mgmi 

Table  3-5.  Comparison  of  the  controllers  (Case  II) 


Controller 

Algorithm 

A;  ^Sl2:S  j  4  ,  c 

robust  statibiliity 

robust  performance 

IW22L 

max  ^(Mji((o,)) 
l<i<t 

ll^ie. 

msm 

~07T8 

1.5154 

0.6209 

Table  3-6.  Comparison  of  the  controllers  (Case  III) 


Controller 

Algorithm 

Aj  =  rdiag[6i,62]:Si.S2  j  A2  c 

robust  statibiliity 

robust  performance 

max  ^(Mji((0,)) 
]<i<( 

l^/L 

max  ^(M(W^)) 
I<iSt 

IML 

Hr 

2.5968 

upper: 

0.3236 

lower; 

0.3236 

3.0904 

upper: 

2.8917 

lower; 

2.7822 

0.3594 

H_ 

1.5154 

upper: 

0.3233 

lower: 

0.3233 

3.0929 

upper; 

1.8896 

lower; 

1.8098 

0.5525 
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Next,  for  Case  II  and  Case  III,  the  block  was  allowed  to  have  two  different  structures 
while  the  full  block  structure  of  ^2  unaltered.  This  can  be  interpreted  as  the  plant  uncertainty 

structure  is  changed  but  the  external  disturbances  and  the  performance  requirements  remains 
unchanged.  In  Case  II,  it  was  assumed  that  A]  ={S!2:S  eCj.  This  means  that  the  uncenain 

parameters  are  only  uniformly  penurbed  diagonally.  In  Case  III  Aj- 

implies  that  the  uncenain  parameters  are  independenUy  perturbed. 


All  three  cases  of  p-analysis  »n  Tables  3-4,  3-5  and  3-6  show  that  both  and 
controllers  possess  reasonably  good  robust  stability  but  not  good  nominal  and  robust  performance. 
However,  the  controller  has  a  better  nominal  and  robust  performance.  On  the  other  hand,  the 
robust  r»erformance  of  the  controller  is  sensitive  to  the  uncertainty  structure  change,  and 
apparently  its  robust  performance  should  improve  if  the  parameters  are  perturbed  uniformly  and 
diagonally. 


The  p-analysis  package  [Balas  et  al.  (1991)]  is  capable  of  calculating  the  upper  and  lower 
bounds  of  p.  As  mentioned  in  Section  2,  the  lower  bound  p  is  not  exact  although  theoretically  an 
equivalence  exists.  The  upper  bound  of  the  exact  p  can  be  reached  for  only  some  specified 
uncertainty  structures  as  seen  in  Table  2-1.  The  upper  bounds  of  the  numerically  computed  p 
match  the  exact  value  of  p  for  AjeC^^^  and  Aj={5l2:SeC}  in  Tables  3-4  and  3-5, 
respectively.  However,  neither  the  upper  bound  nor  the  lower  bound  is  guaranteed  to  match  the 


exact  value  of  p  as  seen  in  Table  3-6  because  for  Aj  = 


Sj 

0 


0 

^2 


■  Sj ,  $2  ^  ^ 


1 

T 


3.3.  p-Synthesis  and  D-K  Iteration 

'fhe.  /fj  optimization  provides  a  unique  controller  while  H_  optimization  uses  an  iterative 
procedure  to  obtain  an  approximate  solution  [Doyle  et  al.  (1989)].  The  p-analysis  could  be  used  to 
guide  the  iterative  process  of  optimization,  and  this  process  is  called  the  D-K  iteration  [Packard 
et  al.  (1993)].  As  described  in  Section  2.5,  any  linear  time-invariant  finite-dimensional  stabilizing 
controller  can  be  used  to  initiate  the  D-K  iteration.  However,  in  the  subsequent  steps  of  D-K 
iteration,  only  synthesis  has  been  used  to  update  the  controller,  and  p-analysis  is  used  to 
evaluate  the  designed  controller  at  each  iteration.  In  the  example  of  this  report,  the  controller  was 
initially  designed  via  H_  optiinization  and,  after  two  passes  of  D-K  iteration,  the  updated  controller 
exhibited  significantly  superior  robustness  properties.  In  each  of  the  two  passes,  a  4th  order 
polynomial  fit  was  chosen  for  the  scaling  matrix  D(s).  The  results  are  summarized  below  in  Tabic 


3-7. 
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Table  3-7.  Comparison  of  the  initial  and  final  controllers  given  by  D-K  iteration 


4y 

42 

D-K 

robust  statibiliity 

robust  performance 

iteration 

number 

I^22L 

max  n(M  1,(03,)) 

l^;L 

max  )i(M((o,)) 
l<i<l 

ll^L 

initial 

0 

1.5154 

I.8IS4  ' 

HH 

1 

1.0175 

wm 

0.6573 

_ 

6.2800 

0.7509 

IHH 

Table  3-7  shows  that  robust  stability  of  the  final  controller  is  significantly  enhanced  since 
is  much  smaller  than  1.  The  nominal  performance  and  robust  performance  are 

also  significantly  improved  since  both  ||iW22lL  i^(M(cd,))  are  made  less  than  1  in  final 

version  of  the  control  system  synthesis. 

The  penurbation  data  obtained  from  the  initial  conirollci  analysis  were  used  to  build  the 
transfer  matrix  Aq{s).  If  the  penurbation  data  obtained  from  the  final  controller  analysis  is  used, 

the  control  loop  closed  by  the  initial  controller  may  become  unstable.  The  rationale  is  that  the  final 
controller  has  better  stability  robustness  than  the  initial  one  so  that  ||d(;(5)|^  of  the  final  controller 

is  greater  than  that  of  the  initial  controller. 

Figures  3-3  and  3-4  show  the  frequency  domain  simulation  of  the  initial  Hoc  controller  and 
final  Hoc  controller,  respectively,  to  compare  their  robust  performance  and  stability.  The  plant 
transfer  matrix  is  G(s)  =G(s)l I  +  '  where  Aq{s)  is  a  stable,  real-rational  transfer 

matrix  constructed  according  to  the  uncertainty  perturbations  obtain'^d  from  the  p-analysis  data. 
Since  a (Aq(s))  ~  uis}))  ^  Vs  c  Sp  where  5/j  is  the  selected  set  o(  frequency  points 
for  numerical  computation  of  the  closed-loop  system  has  the  worst  case  penurbations  among 
the  set  of  actual  plant  transfer  matrices. 

Figures  3-3  and  3-4  compare  the  singular  values  of  the  nominal  loop  transfer  matrices  and 
nominal  sensitivity  matrices,  respectively.  The  loop  with  the  final  controller  has  better 
performance  than  that  with  the  initial  controller  because  the  minimal  singular  values  of  the  loop 
transfer  matrices  of  the  final  controller  is  higher  at  the  low  frequency  range.  Another  interesting 
point  is  that  the  range  of  singular  values  from  maximum  to  minimum  is  much  larger  for  the  final 
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controller  than  that  for  the  initial  controller.  It  is  important  to  note  that  the  conventional  robust 
performance  rule  is  very  conservative  under  the  assumption  that  all  the  singular  values  at  the  low 
frequency  range  must  be  high.  However,  since  the  final  controller  is  constructed  based  on  the  |i 
theory,  it  is  much  less  conservative  than  a  controller  synthesized  via  the  conventional  techniques. 
So  the  robustness  constraints  are  relaxed  and  each  singular  value  can  be  varied  within  a  large  range 
while  the  performance  of  the  final  controller  as  shown  in  Figure  3-3  is  still  higher.  On  the  other 
hand,  the  maximum  singular  value  of  the  loop  transfer  matrix  of  the  final  controller  is  higher  than 
that  of  the  initial  controller.  But  the  robust  stability  of  the  final  controller,  as  shown  in  Table  3-7,  is 
better  than  that  of  the  initial  controller.  This  is  dvie  to  the  same  reason  that  p.-theot7  is  much  less 
conservative  so  that  the  stability  bound  for  the  nominal  plant  can  be  relaxed  while  keeping  the 
stability  robustness  high.  The  example  given  in  this  report  shows  the  difference  between  the 
nominal  stability  and  robust  stability,  and  the  difference  between  the  nominal  performance  and 
robust  performance.  For  multivariable  system  controller  analysis  and  synthesis,  we  should  not 
only  consider  the  nominal  stability  and  performance  but  also  the  robust  stability  and  performance. 
However,  the  benefit  of  the  robustness  is  gained  at  the  expense  of  an  increased  order  of  the 
controller  in  the  D-K  iteration.  For  real  time  implementation,  a  reduction  of  the  control  lei  model 
order  (Balas  et  al.  (1991)]  might  be  necessary. 
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Figure  3.3  Loop  Transfer  Matrices  of  the  Control  System 


Singular  Value  Of  The  Loop  Sensitivity  Matrix 


Figure  3.4  Loop  Sensitivity  Matrices  of  the  Conirol  System 
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3.4.  The  Control  Analysis  And  Synthesis  Procedure  Using  ^ 

Procedures  for  numerical  computations  in  |i-analysis  and  p-synthesis  are  discussed  in  view 
of  the  analytical  results  of  Section  2  and  the  simulation  results  of  Section  3.  Given  a  nominal 
plant,  the  perturbed  parameter  and  uncertainty  structure,  the  design  specifications,  and  a  set  of 
MIMO  controllers  (that  are  already  designed),  the  following  procedural  rules  are  presented  for 
analysis  and  comparison  of  these  controllers  in  view  of  their  nominal  and  robust  stability  and 
performance. 

Step  1:  Break  the  loop  at  selected  points  according  to  the  uncertainty  structure  and  perturbed 
parameters; 

Step  2:  Select  the  weighting  functions  for  stability  and  performance  according  to  the  design 
specifications  and  knowledge  of  the  uncertainty  structures; 

Step  3:  Build  the  augmented  open-loop  system  structure  which  includes  the  nominal  plant,  and  the 
weighting  functions; 

Step  4;  Use  each  of  the  designed  controllers  to  close  the  loop  and  construct  the  augmented  M-block 
for  p-analysis; 

Step  5;  Check  whether  each  pole  of  the  M-bkKk  is  located  in  the  left  half  s-plane.  Choose  those 
controllers  that  arc  stable. 

Step  6;  Use  the  p-analysis  program  (MU_ANLY.M  described  in  Section  4)  to  compare  ||W22li^ 

for  each  design  and  choose  those  controllers  that  most  closely  satisfy  the  nominal 
performance  requirements. 

Step  7:  Compare  )llo„  for  each  design  and  choose  those  controllers  that  most  closely  satisfy 

the  robust  stability. 

Step  8;  If  the  robust  performance  is  the  major  concern  of  the  design  objective,  then  )|1^  of  the 

candidate  controllers  should  be  compared  first.  Otherwise,  use  different  uncertain.^ 
structures  and  repeat  p-analysis  to  verify  if  the  robust  stability  of  the  controller  is  sensitive 
to  a  change  in  the  uncertainty  structure.  Choose  those  conmollers  that  most  closely  satisfy 
the  robust  performance  or  arc  relatively  less  sensitive  to  changes  in  the  uncertainty- 
structure. 

Step  9:  Conduct  a  time  domain  simulation  on  the  controllers  that  have  passed  the  above  tests  for  the 
actual  closed-loop  system  structure  under  the  worst  case  perturbations.  Select  those 
controllers  which  provide  a  good  performance  of  the  closed-loop  system. 
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4.  DESCRIPTION  OF  THE  COMPUTER  CODE 

A  number  of  programs  are  written  in  the  MATLAB  environment  [Moler  et  al,  (1987);  Balas 
et  al.  (1991)]  for  robustness  analysis  with  any  multivariable  plant  model  and  uncertainty  blocks. 
The  program  MU_EXAMP.M  demonstrates  the  usage  of  the  above  programs  by  analyzing  the 
flight  control  system  of  an  advanced  aircraft,  discussed  in  Chapter  3  and  following  the  program 
structure  shown  in  Figure  4- 1 .  However,  to  use  these  programs  in  the  specific  p-analysis  and  |i- 
synthesis  problem,  the  user  sliould  write  his  or  her  own  program. 

MUINIT.M 

This  program  creates  an  interconnection  stmcture  of  a  given  open  loop  system.  It  must  be  run 
before  any  controller  design  and  p-analysis. 

Input:  A,  B,  C.  D  -  nominal  plant  state  space  matrices. 

Gdel  ■  gain  of  uncertainty  weighting  function. 

Gp  -  gain  of  performance  weighting  function. 

Output:  Wdel  -  uncertainty  weighting  function. 

\Vp  -  pcrfomiancc  weighting  function. 

p_block  -  interconnected  strucnire  of  the  open  loop  system. 

Since  the  interconnection  structure  is  specific  for  each  individual  system,  the  program  written 
inside  MU_INiT,M  should  be  changed  in  each  case.  The  user  should  know  how  to  build  the 
interconnection  structures  using  the  MATLAB  function  "sysic"  [Balas  et  al.  (1991)].  Figure  4-2 
gives  the  interconnection  structure  used  in  the  example  given  in  Chapter  3.  Three  different 
methods,  namely,  LQG,  H2  and  Hoo  for  the  initial  step  of  controller  design  are  provided  in  the 
programs  called  MU_L(3G  M.  MU_II2.M  and  MIJ_HINF.M,  icspectively. 
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MU  LQG.M 

This  program  carries  out  LQG  controller  design.  It  can  be  used  as  the  initial  step  of  the  p 
synthesis. 

Input;  A,  B,  C,  D  -  plant  state  space  matrices 

Qc,  Rc  -  weighting  matrices  for  LQR  design 

Qf,  Rf  -  plant  disturbance  and  measurement  error  covariance,  respectively 
Output:  F  -  controller  gain 
L  -  observer  gain 
k_comp  -  compensator  structure 

M_plant  -  closed-loop  system  structure  as  shown  in  Figure  4-2, 


Figure  4-2.  The  Closed-Loop  System  M_plant. 

MU_H2.M 

This  program  carries  out  H2  controller  design.  It  can  be  used  as  the  initial  step  for  p 
synthesis. 

Input:  p_block  -  open-loop  system  structure 

Output:  k_comp  -  compensator  structure 

M_plant  -  closed-loop  system  structure 
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MU  HINF.M 


This  program  does  Hoo  controller  design.  It  can  be  used  as  the  initial  step  for  |i-synihesis. 

Input:  p_block  -  open-loop  system  structure 

Output;  k_comp  -  compensator  structure 

M_plant  -  closed-loop  system  structure 

The  next  step  after  controller  design  and  evaluation  is  |i-analysis  as  seen  in  the  flow  chart  of 
Figure  4-1. 


Figure  4-3.  The  Structure  Used  in  the  p-3nalysis 

MU  ANLY.M 


This  program  carries  out  |i-analysis. 

Input:  M_plant  -  closed-loop  system  structure 

blkl  -  block  structure  of  Ai  as  shown  in  Figure  4-3. 
blk2  -  block  structure  of  A2  as  shown  in  Figure  4-3. 


blk  -  block  structure  of  A 


A,  0 
0  A2 


Output;  frequency  domain  plot  of  and  ||M22!1^  as  the  measure  of  nominal  performance. 


frequency  domain  plot  of  max(/i(Mj  j)),  and  |jAj|l^  as  the  measure  of  robust 

stability 

frequency  domain  plot  of  max(/i(M)),  and  ||A||^  as  the  measure  of  robust 

performance. 


At  this  stage  the  user  may  choose  to  further  improve  the  controller  via  the  D-K  algorithm  of  |i- 
synthesis  or  directly  examine  the  performance  of  controller  via  time-domain  simulation. 
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Note:  When  typing  the  input  of  blkl,  blk2,  and  blk  from  the  keyboard,  it  is  required  to  use  the 

same  format  as  prescribed  in  [  Balas  et  al.  (1991)  pp.  4-70]. 


MU  DK.M 


This  program  dose  D-K  iteration  as  a  part  of  |i-synthesis. 

Input:  p_block  -  original  open-loop  system  structure 

rpdvec  -  D  matrices  obtained  by  |i-analysis 
rpsens  -  sensitivity  function  of  p.(*)  obtained  by  |i-analysis 


blk  -  block  structure  of  4  = 


0 


0 

Al 


Output:  k_comp  -  compensator  structure  with  increased  order  during  the  iteration.  This  is 
designed  according  to  the  new  open  loop  system  structure  shown  in  Figure  4-4. 

M_plant  -  closed-loop  system  structure  with  increased  order  during  the  iteration 


Figure  4-4.  Augmention  via  D-Matrices  in  the  D-K  Iteration 
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MU  SIMU.M 


This  program  does  time  domain  simulation  of  the  closed  loop  system  with  the  worst  perturbation 
as  a  stable  transfer  matrix. 

Input:  M_plantl  -  closed-loop  system  structure  of  the  initial  design 

in_pvec  -  perturbations  of  the  initial  design  obtained  from  |a,-analysis 

in_bnds  -  upper  and  lower  ji  bounds  of  the  initial  design  obtained  from  u-analysis 

M_plant2  -  closed-loop  system  structure  of  the  final  design 

fin_pvec  -  penurbations  of  the  final  design  obtained  from  p-analysis 

fin_bnds  -  upper  and  lower  p.  bounds  of  the  final  design  obtained  from  p-analysis 

n_block  -  number  of  blocks  to  be  approximated  by  transfer  matrix 


blk  -  block  structure  of  A  = 


0 


“I 


Output:  M^pertl  -  closed-loop  system  structure  of  the  initial  design  with  the  worst  perturbation  as 
shown  in  Figure  4-5. 

M_pert2  -  closed-loop  system  structure  of  the  final  design  with  the  worst  perturbation 
u  -  time  domain  input  signal 
yl  -  output  signal  of  the  initial  design 
y2  -  output  signal  of  the  final  design 


M_pen(s) 


Figure  4-5.  Closed-Loop  System  with  the  Worst-Case  Perturbation 


29 


NAWCADWAR-93010-60 


MU  FRQ.M 

This  program  plots  singular  values  of  loop  transfer  matrices  and  sensitivity  matrices  for  a 
closed-loop  system. 

Input:  ap,bp,cp,dp  -  plant  state  space  matrices 

ak,bk,ck,dk  -  compensator  state  space  matrices 
ak2,bk2,ck2,dk2  -  second  compensator  state  space  matrices 

n_comp  -  number  of  closed-loop  systems.  If  n_comp=l,  then  only  one  loop  is  plotted;  if 
n_comp  =2,  then  two  loop  functions  ai-e  plotted  in  one  graph. 

Output:  sxl  -  singular  values  of  the  loop  transfer  matrices 

syl  -  singular  values  of  the  loop  sensitivity  matrices 

sx2  -  singular  values  of  the  second  loop  transfer  matrices,  if  n_comp=2 

syl  -  singular  values  of  the  loop  sensitivity  matrices,  if  n_comp=2 
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5.  SUMMARY  AND  CONCLUSIONS 

This  report  presents  the  results  of  control  systems  analysis  and  their  applications  to  control 
systems  synthesis  for  high  performance  aircraft,  and  provides  a  proven,  tractable  and  portable 
methodology  to  quantify  the  stability  margins  and  robustness  of  highly  integrated  digital  fly-by¬ 
wire  control  systems.  The  Hoo-based  ji-analysis  approach  has  been  recommended  for 
quantitative  evaluation  of  robust  stability  and  performance  of  multi-input  multi-output  (MIMO) 
control  systems.  This  approach  for  control  systems  evaluation  offers  an  alternative  to  the  notions 
of  conventional  gain  margin  and  phase  margin. 

A  survey  of  recent  literature  has  shown  three  major  different  approaches  to  quantitatively 
evaluate  the  robust  stability  and  performance  of  MIMO  control  systems.  The  first  approach  relies 
on  the  structured  singular  value  (|i)  and  Hoo-analyses  in  the  framework  of  Linear  Fractional 
Transformations  (Lf’O  of  MIMO  systems  that  are  subjected  to  uncertainties  with  known  bounds. 
The  sero.-i'i  approach  generalizes  the  definition  and  calculation  of  gain  and  phase  margins  for 
MIMCj  -'^y.su.ms  analogous  to  those  for  single-input,  single-output  (SISO)  systems.  The  third 
approac'  examines  the  relationships  between  the  open  loop  gain  and  phase  with  respect  to  closed 
.oc^  jitivity  and  complementary  sensitivity  functions  for  SISO  systems,  and  then  formulates 
heuiistic  rules  for  shaping  the  loop  transfer  functions  for  MIMO  systems.  The  first  approach  of  1.1- 
analysis  for  robustness  assessment  of  MIMO  systems  constitutes  the  main  body  of  this  report 
because  it  has  significant  advantages  over  the  remaining  two  approaches  as  discussed  below. 

The  second  approach  attemptr  to  generalize  the  conventional  concepts  and  computation 
algorithms  of  gain  margin  and  phase  margin  for  MIMO  systems.  The  concepts  of  gain  margin  and 
phase  margin  were  initially  based  on  gain  parameter  perturbations  of  a  given  controller  for  SISO 
closed-loop  systems.  The  algorithm,  generalized  for  MIMO  systems  in  the  second  approach,  can 
only  deal  with  parameter  variations  of  a  constant  gain  matrix.  Since  a  MIMO  plant  is  generally 
represented  by  a  stabilizable,  real-rational  transfer  matrix,  the  system  gain  margin  is  not  uniquely 
defined.  Therefore,  analysis  of  MIMO  control  systems  based  on  gain  margin  and  phase  margin  is 
not  a  straight- forw'ard  extrapolation  of  that  of  SISO  systems.  In  contrast,  the  concept  of  structured 
singular  value  (p)  deals  with  all  types  of  bounded  parametric  and  non-parametiic  uncertainties. 
Apparently  p-analysis  appears  to  be  much  more  general  than  any  of  the  reponed  algorithms  for 
gain  margin  and  phase  margin  evaluation.  The  p-analysis  can  handle  not  only  the  parameter 
variations  in  the  dynamic  controller  K(s)  but  also  the  uncertainties  in  the  plant  model  and 
parameter  perturbations  anywhere  in  the  loop.  To  summarize,  p-analysis  is  capable  of  quantifying 
the  stability  and  robustness  of  the  closed-loop  system  under  any  destabilizing  factor.  Funhermore, 
the  algorithms  in  the  second  approach  are  relatively  more  conservative  because  singular  values  of 
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the  transfer  matrices  are  often  used  as  the  norm  to  bound  the  perturbations,  and  the  parameters  are 
only  allowed  to  vary  diagonally.  In  contrast,  p-analysis  does  not  suffer  from  these  problems 
because  the  notion  of  structured  singular  value  is  directly  based  on  the  Nyquist  stability  theorem. 

Although  the  third  approach  which  is  based  or  loop  transfer  matrices  is  more  general  than  the 
second  approach  for  MIMO  control  systems,  it  only  provides  a  qualitative  description  of  the 
stability  and  performance  of  closed-loop  systems.  Apparently,  quantitative  analysis  of  stability  and 
perfon  nee  can  be  reliably  carried  out  by  the  algorithms  based  on  p.  Therefore,  p-analysis  is 
conside  ?d  to  be  the  most  effective  tool  for  control  systems  evaluation  as  reported  in  open 
literature 

Results  of  simulation  experiments  using  the  flight  dynamic  model  of  an  advanced  aircraft 
have  been  presented  to  demonstrate  how  the  nominal  and  robust  stability  and  performance  can  be 
evaluated  via  p-analysis.  Furthermore,  it  is  shown  that  both  stability  and  robustness  can  be 
improved  using  the  tools  of  p-synthesis.  A  software  package  has  been  developed  in  the 
environment  of  MATLAB  for  analyzing  robustness  of  MIMO  control  systems.  This  package  is 
built  upon  the  commercially  available  p- Analysis  and  Synthesis  Toolbox  (Balas  et  al.  (1991)]. 
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APPENDIX  A 

The  Principle  of  Multivariable  Nyqulst  Stability  Criterion 

Let  N(A,f(s),C)  denote  the  number  of  clockwise  encirclements  of  the  point  A  in  the 
complex  plane  by  the  image  of  the  clockwise  contour  C  under  the  mapping Using  this 
notation,  the  principle  of  the  argument  can  be  written  as: 

N(0,f(s).C)  =  Z-P 

w'here  Z  is  the  number  of  zeros  of f(s)  inside  C,  and  P  is  the  number  of  poles  of f(s)  inside  C. 

Suppose  f(s)  =  f\(s)f2(s)  and  there  are  no  pole-zero  cancellations  between  fi(s)  and 
f2(s)  Then,  the  Nyquist  stability  criterion  is  stated  as  follows: 

N(Q.fi  \s)f2  (s).  C)  =  N(0.  /i  (s).  C)  +  NrO,  /jf  5).  C) 

Proof  of  the  Nvauist  Stability  Criterion: 

Let  the  contour  Df^  enclose  right  half  plane.  The  closed  loop  system  is  stable  if  and  only  if 
N(0,<i>,i(s).Df()  =  Q 

Therefore,  N(Q,<t>(,i(s)det(l  ■¥G(s)),D/i)  -  N(0,<t>pi(s),D/^)  +  N(0,det(I  +  G(s)),Dfi)  =o. 

Let  P^^  gi  denote  the  number  of  unstable  open-loop  poles  of  G(s).  Then 
N(0.<l>^l(s}.Dii)=P^^^l  and  N(0,det(I -^-Gis^.D^)  = -P^  ,,i 


Figure  A-1 .  Multivariable  Feedback  System 

Given  the  plant  in  Figure  A-1,  the  relationship  between  the  open  loop  and  closed  loop 
characteristic  polynomials  is  as  follows  [Desoer  ct  al.  (1977)]; 
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<Pd  =<^oldetll  +  G(s)] 

(A-1) 

or 

0d  =  del  Is!  -AdJ 

(A-2) 

where  =  detlsl  -  A]  is  the  open  loop  characteristic  polynomial,  A  and  Ad  are  the  system 
matrices  of  the  open  loop  system,  G(s),  and  the  closed  loop  system,  G(s)[I  +G(s)]~^ , 
respectively. 


The  Nyquist  stability  criterion  states  that  the  closed-loop  system  is  stable  if  and  only  if  the 
number  of  counter-clockwise  encirclements  of  the  origin  by  the  image  of  the  clockwise  Nyquist 
contour  under  the  mapping  det(l  -t-  G(s))  equals  the  number  of  unstable  open-loop  poles  of  the 

open-loop  transfer  matrix  G(s).  The  proof  of  this  criterion  is  based  on  the  theory  of  complex 
analysis  as  outlined  earlier  in  this  Appendix.  Several  researchers  [Yeh  et  al.  (1985a),  Molander  and 
Willems  (1980),  Tao  et  al.  (1991),  Lehtomaki  et  al.  (1981)]  have  attempted  to  extend  the  concepts 
of  gain  and  phase  margins  of  SISO  systems  to  address  the  stability  robustness  of  MIMO  systems. 
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APPENDIX  B 

APPROACH  2:  Multivariable  Gain  Margin  and  Phase  Margin 

Now  we  introduce  the  pertinent  definitions  and  some  of  the  useful  algorithms  for  the  afore¬ 
mentioned  two  approaches  of  robustness  assessment. 

Some  diversity  and  ambiguity  exist  in  defining  the  MIMO  gain  and  phase  margins.  The 
common  practice  of  one-loop-at-a-time  stability  margins  fails  to  account  for  the  simultaneous 
variations  in  MIMO  feedback  systems  and  hence  may  not  be  acceptable  as  relative  stability 
measures.  The  norm-bounded  robustness  criteria  guarantee  closed-loop  stability  but  they  do  not 
explain  how  the  individual  elements  of  the  control  gain  matrix  may  be  varied  without  destabilizing 
the  closed-loop  system.  To  this  effect,  Yeh  et  al  (1985a)  introduced  two  different  definitions  of 
stability  margins. 

Definition  B-1:  independent  gain  margins  are  limits  within  which  the  gains  of  all  feedback  loops 
may  be  varied  simultaneously  and  independently  without  destabilizing  the  system  whib  the  phase 
angles  remain  at  their  nominal  values.  Independent  phase  margins  are  limits  within  which  the  phase 
angles  of  all  feedback  loops  may  be  varied  simultaneously  and  independently  without  destabilizing 
the  system  while  the  gains  of  all  feedback  loops  remain  at  their  nominal  values.  However,  these 
stability  margins,  if  computed  via  the  singular-value  based  robust  stability  criteria,  tend  to  be  very 
conservative. 

Definition  B-2:  Uniform  gain  margins  are  limits  within  which  the  gains  of  all  feedback  loops 
may  be  varied  uniformly  without  destabilizing  the  feedback  system  while  the  phase  angles  remain 
at  their  nominal  values.  Uniform  phase  margins  are  limits  within  which  the  phase  angles  of  all 
feedback  loops  may  be  varied  uniformly  without  de.stabilizing  the  system  while  the  gains  remain  at 
their  nominal  values. 

Remark::  A  major  difference  between  MIMO  and  SISO  systems  is  that  the  stability  margins  of 
MIMO  systems  vary  relative  to  the  point  where  the  complex  loop  gains  are  measured.  For  a  general 
feedback  system,  if  the  gain  perturbations  arc  calculated  at  input  of  the  plant,  then  the  simultaneous 
perturbation  in  each  loop  may  be  represented  by  a  diagonal  perturbation  matrix  L(s)  preceding  the 
plant  G(s). 

For  independent  margins:  L{joi)  -  cliagl hj( a})e^‘^ h2( . (B-1) 

For  uniform  margins:  H  jO))  =  [i( (B-2) 
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There  are  several  different  ways  to  calculate  the  stability  margins  (Ych  et  al.  (1985a), 
Molander  and  Willems  (1980),  Tao  et  al.  (1991)].  However,  we  discuss  the  following  three 
algorithms  in  this  report. 

Robustness  Bound  For  Penurbaiion  Errors:  [Yeh  et  al.  (198ja)) 

Let  G(s)  and  G(5)  represent  the  nominal  and  perturbed  open-loop  transfer  matrices, 
respectively,  and  the  multiplicative  and  inverse  multiplicative  perturbations  are  given  as  follows: 


[  [G(s) -G(s)]g~  (s)  for  multiplicative  perturbations 

t(s)  =  'i  .1  _  J 

[[G  (s)  -  G  (s)]G(s)  for  inverse  multiplicative  perturbations 


(B-3) 


In  eqs.  (B-1)  and  (B-2),  the  gain  and  phase  perturbations  L(s)  =  E(s)  +  I.  The  singular  value 
based  robusmess  tests  are  given  as  follows  (Lehtomaki  et  al.  (1981)]: 


a[L  /]<a^/-(-G  'j  for  multiplicative  perturbation,  and 

-  /j  <  a  <  a[/  +  G],  a  ^  I  for  inverse  multiplicative  perturbation. 


(B-4) 

lB-5) 


It  has  been  shown  [Yeh  et  al.  (1985b)]  by  using  the  inverse  Nyquisi  criterion  and  the  Nyquist 
criterion  that  the  robustness  tests  for  (B-3)  are: 


max  min  r  .-ii 

,  |AaL-/J(<  A,|/  +  G  I  for  multiplicative  perturbation,  and  (B-6) 

^^WLX|,  -7  I  min.  , 

.  |A,yL  ^  -  /y|  <  a  <  .  \X^[!  +  GjJ.  0  <  a  <  ; 


for  inverse  multiplicative  perturbation. 


(B-7) 


The  independent  gain  margin  (IGM)  is  obtained  by  substituting  (B-2)  into  (B-4)  and  then  letting 
the  phase  angles  9]  =  62  =  •=  9^  =  0,  and  ai  + 


1-aj  <  Pi(0})  <  l-t-aj 


Vi 


(B-8) 


Letting  the  gains  =  P2  Pn~^<  independent  phase  margin  (IPM)  is  obtained  as: 

-25in~‘(fli/2)  <  0,(0))  <  2jin"*(aj/2)  Vi,  a\<l  (g 

and  ~n<9p.a))<n  Vi,  a,  >  2  (p.lO) 


36 


NAWCADWAR-93010-60 


Similarly,  substituting  (B-3)  into  (B-5)  and  the  !GM  and  IPM  for  inverse  multiplicative 
perturbations  are  obtained  as: 

— ^ — <P:{(0)< — - —  Vj  where  aj  =  of/ +  GJ,  and 

1  +  a,  '  l-cti  ^  ^  (B-11) 

and  -2sin~^{a]l2)  <  6,(03)  <2sin~^{ail2)  Vi  (B-12) 


Substituting  (B-3)  into  (B-6)  and  (B-7),  the  uniform  gain  margin  (UGM)  and  uniform  phase 
margin  (UPM)  can  be  obtained  as: 


For  multiplicative  penurbation:  Oq  = 


min 

i 


under  the  following  constraints: 


1  -oq  <  p(a))  <1  +  ^:0 

-2sin~^{aQl2.)  <  6(0))  <  2sin~^{aQ/2)  % 

-n  <  6(03)  <  n  olq>2 


(B-13) 

(B-14) 

(B-15) 


For  inverse  multiplicative  penurbation:  «o  ” 


min 

i 


|A,(;  +  G]j 


under  the  following  constraints: 


- <  p(0))  < - 

l-t-ttQ  1-ao  -  (B-16) 

-25in~\aQ/2)  <  6(03)<  2sin'\aQl2) 


The  advantages  of  using  this  algorithm  to  calculate  M!MO  gain  margins  and  phase  margins  are  that 
lelatively  simple  formulae  can  be  derived  from  the  robustness  bound  concepts.  Although  the 
independent  gain  and  phase  margins  are  rather  conservatively  based  on  singular  values,  the 
combined  use  of  IGM,  IPM,  UGM  and  UPM  yields  a  less  conservative  estimate  of  the  stability 
perturbing  region  in  the  gain  space  and  phase  space  {Yeh  et  al.  (1984)].  How-ever,  the  major 
disadvantage  is  that  these  perturbations  are  only  constrained  at  diagonal  elements  of  the 
perturbation  matrix  There  is  apparently  no  discussion  about  the  off-diagonal  tenms. 

Perturbation  With  Nonlinear  Cone  Boundary  (Molander  and  Willems  (1980)] 

For  MIMO  linear-quadratic  regulators,  there  is  Kalman's  inequality  which  implies  that  this 
kind  of  controllers  possess  a  60°  phase  margin,  infinite  gain  margin,  and  50  percent  gain  reduction 
tolerance  [Safonov  and  Athans  (1977)].  Using  this  result,  Molander  and  Willems  (1980) 
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developed  a  procedure  for  synthesizing  state  feedback  control  of  multivariable  systems  with 
specified  gain  and  phase  margins.  Considering  the  linear  time-invariant  finite-dimensional  systems: 

x-Ax^Bu  (B-18) 

with  xeR'' ,  ueR"^,  and  the  pair  (A,  B)  being  controllable,  and  the  desired  input  being  given  by 
the  linear  time-invariant  state  feedback  control  law: 

u  =  -llx  (B-19) 

The  above  equation  (B-19)  is  called  a  robust  control  law  with  robustness  sector  [K],K2)  with 

->»  <  7  <  if  the  closed-loop  system  is  globally  asymptotically  stable.  The  sector 

condition  is  illustrated  in  Figure  B-1  for  the  single-input  case.  The  control  input  is  given  by 

u  =  -f(Jx.t)  (B-20) 


with  f:R  xR—^R  being  any  nonlinear  function  satisfying  the  following  condition: 


sup 

<T,t 


[f(o.i)  -  Kic^f  (/( cT.rj  -  K20)  <  0 


sup\f(a.t)-\l2{Kx^K2)<^  1,^ 

Or,  equivalently, - p- - 1<— 

2  ' 


o,t 


(B-21) 


(B-22) 


Figure  B-1,  An  Illustration  of  the  Sector  Condition 
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The  robustness  requirements  also  allow  for  an  interpretation  in  terms  of  gain  and  phase 
margins  and  gain  reduction  tolerance  (defined  below).  The  control  law  (B-19)  has  a  gain  margin  « 
ig>\)  if  the  closed-loop  system  is  globally  asymptotically  stable  for  all  control  laws: 

u  =  -All  X  (B-73) 

where  A  is  any  symmetric  matrix  with  eigenvalues  1  <  A,  <  g .  The  system  is  said  to  have  a  phase 
margin  (0  <  ^  <  ;r)  if  it  is  globally  asymptotically  stable  for  all  control  laws  (B-23)  with  A 
being  any  unitary  matrix  possessing  eigenvalues  A,  =  exp{j(p^)  with  |0il  ^  0.  It  is  said  to  have  a 

gain  reduction  tolerance  p  if  it  is  globally  asymptotically  stable  for  all  control  laws  (B-23)  with  A 
any  symmetric  matrix  with  eigenvalues  (1  -p)  <  A,  <  /.  The  relationship  between  robustness 
sector  {K],K2)  and  the  gain  and  phase  margins  are: 


8  =  ^2 

(B-24) 

p  =  (\-K^)xlOO 

(B-25) 

COS(  0)  =  — — = - 

K^  +  K2 

(B-26) 

Using  the  multivariable  control  system  stability  theory,  a  methodology  has  been  developed 
[Molander  and  Willems  (1980)]  to  synthesize  the  linear  quadratic  control  law  under  prescribed  gain 
and  phase  margins.  A  similar  synthesis  procedure  has  been  developed  for  the  discrete-time  state 
feedback  control  law  by  [Lee  and  Lee  (1986)]. 

The  advantages  of  the  above  control  algorithm  are  that  the  synthesis  is  based  on  the 
framework  of  the  general  state  feedback  control  law,  which  can  be  directly  obtained  with  a  given 
robustness  sector  by  solving  the  matrix  Riccati  equation.  'Fhe  major  disadvantage  of  this  algorithm 
is  that  it  is  apparently  restricted  to  the  full-state  feedback  control  case  and  does  not  address  the 
general  compensator  structures. 

Unitv-Feedback  Multivariable  Control  System  ITaoet  al.  (1991)] 

In  this  method,  the  boundaries  of  gain  and  phase  margins  of  multivariable  control  systems 
are  obtained  by  decomposition  into  augmented  open-loop  transfer  functions.  Then  the  characteristic 
equations  are  formulated  and  the  stability  margins  are  tested.  Consider  the  multivariable  feedback 
control  systems  shown  in  Figure  B-2. 
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R(':Hr,(s)r  (s)  ...r(s)]  ^ 

I  2  m 


Y(s)=[y  (s)y(s) ...  y(s)]^ 
12  m 


I 


Figure  B-2.  Multivariable  System  Represented  by  Augmented  Open-Loop  Transfer  Function 


The  proper  transfer  function  G(s)  is: 

'Gn 

Gn  • 

•  Gu  . 

G\m 

Gii 

G22 

G2i 

G2m 

G(s)  = 

Gn 

Gn  • 

■  Gu  . 

Gim 

=  N(s)D~\s) 

Gm\ 

Gml  • 

■  G^i  . 

G/nn  j 

where  the  relatively  coprime  polynomial  matrices, /V('5>  and  D(s),  are  er^pressed  as: 
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"12 

••  "1( 

"Im 

"21 

"22 

"2i 

"2m 

N(s)  = 

"jl 

"i2 

••  "im 

."ml 

"mi 

•••  "mn 

^c\ 

0 

0  .. 

0  ■ 

0 

^c2 

0 

0 

D(s}  = 

0 

0 

• 

0 

0 

0 

.  0  .. 

^cn_ 

(B-28) 


(B-29) 


The  characteristic  equation  of  the  closed  loop  system  is 


dc\  +  ni  1 

"12 

nu 

"im 

"21 

^c2  +  "22 

nil 

nim 

det[D(s)  +  N(s)]  =  del 

"il 

"i2 

..  dci+na  . 

nim 

1 

"ml 

"m2 

tlfni 

(B-30) 


Applying  the  determinant  property  (Kreyszig  (1972)),  eq.  (B-30)  can  be  decomposed  into  the 
following  m  different  expressions 

det{D(s)-^  N(s)]  =  der[Do,(5)]  +  der[Ao,(5))  =  0,  /  =  l,2,...,m  {B-31) 

where  dei[NQi(s)]  det[DQi(s)]  are  defined  as; 


dc\  +"11 

"12 

..  n,,  .. 

"im 

"21 

^c2  +  "22 

nil 

nim 

=  det 

"il 

"i2 

..  na 

nim 

"ml 

"m2 

n/ni 

■■  dcn^n, 

=  0 


(B-32) 
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dct[DQ^(s)]  =  de! 


+  ”1 1 

"12 

..  0  . 

«21 

^c2  +  ”22 

0 

”2m 

-  4.  • 

^im 

”ffi2 

..  0  . 

=  0 


(B-33) 


Now  the  augmented  open-loop  transfer  function  (AOLTF)  is  defined  as: 


80i(‘^)  = 


_  der[%('^)]  ^  Ni(s) 
dei[DQ^(s)]  D,(s) 


(B-34) 


For  a  unity-feedback  mxm  multivariable  system,  there  are  m  AOLTFs.  Each  of  the  goi(s)  is  a 
transfer  function  of  yi(s)^r^(s)  effected  by  the  ith  (open)  unity-feedback  loop,  and  can  be 
considered  as  an  open-loop  transfer  function  of  an  unit  feedback  control  system  shown  in  Figure 
B-2.  The  gain  and  phase  margins  of  an  unity-feedback  multivariable  control  system  are,  in  fact,  the 
gain  margins  (GM)  and  phase  margins  (PM)  of  all  AOLTFs.  For  an  mxm  unity-feedback 
multivariable  control  system,  there  arc  m  sets  of  gain  margins  and  phase  margins.  The  calculation 
of  each  set  is  based  upon  the  corresponding  AOLTF  and  can  be  checked  by  Nyquist  plots. 

The  advantage  of  the  above  synthesis  algorithm  is  that  the  plotting  methods  used  in  SISO 
feedback  loop  stability  tests  can  be  used  directly  and  therefore  the  gain  and  phase  margins  of  each 
AOLTF  are  easily  formulated.  The  relationship  of  tlie  stability  margins  of  AOLTF  with  parameter 
perturbations  of  the  original  system  is  discussed  by  Tao,  Cheng  and  Han  (1991).  These  parameters 
could  be  the  original  system  gain  perturbations.  However,  funher  work  is  needed  is  to  find  the 
direct  relationships  of  the  stability  m,,rgins  of  AOLTFs  with  the  original  control  system. 

Each  of  the  above  three  algorithms  is  airr.“d  at  calculation  of  the  MIMO  gain  and  phase 
margins  analogous  to  those  of  SISO  systems.  However,  the  ultimate  goal  is  to  shape  the  closed- 
loop  functions  such  that  it  satisfies  the  performance  and  stability  robustness  requirements.  Since 
there  is  a  direct  relationship  between  the  open  loop  gain  and  phase  and  closed-loop  characteristics 
for  SISO  systems,  the  gain  margin  and  phase  margin  are  two  parameters  to  be  adjusted  to  shape 
the  closed-loop  transfer  function.  There  could  be  other  parameters  to  adjust  the  closed  loop  transfer 
functions  for  MIMO  systems.  Following  this  concept,  Freudenburg  and  Looze  (1988)  have 
investigated  the  MIMO  closed-loop  system  properties,  and  developed  rules  for  shaping  the  closed- 
U  )p  sensitivity  function  and  complementary  sensitivity  function  by  adjusting  the  open  loop  gain, 
phase  and  directions. 
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APPENDIX  C 

APPROACH  3:  Multivariable  Loop  Shaping 

Classical  "loop-shaping"  design  methods  proceed  directly  with  manipulation  of  open  loop 
gain  and  phase  to  alter  the  feedback  propenies  of  the  system  [Freudenburg  and  Looze  (1986), 
Freudenburg  (1990),  Freudenburg  and  Looze  (1988)].  Clearly,  one  reason  for  the  success  of 
these  methods  is  that,  for  a  SISO  system,  open  loop  gain  and  phase  can  be  readily  related  to  the 
feedback  properties.  The  loop  transfer  function  L( JO))  has  the  following  properties: 

\S(jQ})\  «  1 

\L(jci))\»  I  and 

T(ja))  =  l  (C-1) 

and 

\T(jfi))\  «  I 
\L(jo))\  «  1  <=>  and 

S(jO))  -  1  (C-2) 

At  frequencies  for  which  the  open  loop  gain  is  approximately  unity,  feedback  propenies  depend 
critically  upon  the  open  loop  phase: 

|L(ycu)|*l  |5(ycoj|»l 

and  <=>  and 

zL(j(o)^±m^  lr(yco)i»i  (C-3) 

Heuristic  rules  used  in  design:  First,  large  loop  gain  >i,,.lds  small  sensitivity  and  good  disturbance 
rejection  properties,  although  noise  appetirs  directly  in  the  system  output.  Second  small  loop  gain 
is  required  for  small  noise  response  and  for  robustness  against  large  multiplicative  uncertainty. 
Finally,  at  frequencies  near  gain  crossover  (i.e.,  \L(jUi)\  =  1),  phase  of  the  open  loop  system  must 

remain  bounded  sufficiently  far  away  from  ±180^  to  provide  an  adequate  stability  margin  and  to 
prevent  amplifying  the  disturbances  and  noise. 

Generalization  of  Scalar  Concepts  to  Multivariable  Systems 

Under  certain  conditions,  multivariable  systems  behave  sufficiently  like  SISO  systems  such 
that  classical  results  can  be  easily  extended.  First,  it  is  assumed  that  the  plant  to  be  controlled  has  a 
well-conditioned  transfer  function  matrix.  It  is  also  necessary  to  assume  that  the  disturbances  are 
injected  into  the  system  at  the  same  loop-breaking  point  at  which  uncertainties  exist,  and  that  the 


43 


NA\VCAD\VAR-93010-60 


levels  of  uncertainty,  noise,  and  disturbances  are  approximately  equal  in  iill  loops  of  the  system. 
Hence,  die  disturbance  and  noise  signals  should  satisfy  the  following  conditions  of  upper  bound; 

\d(jco}\<  Mdico)  Vco  and 

\n(i03)\<  M,j{co)  V(o 

where  Mjoj)  is  an  upper  bound  of  the  norm  of  a  given  transfer  matrix  x(j(jo). 

We  assume  that;  (i)  each  of  the  outputs  of  the  system  is  significant;  and  (ii)  the  response  of 
each  to  disturbances  and  noise  is  desired  to  be  small.  Therefore,  the  following  assumptions  need  to 


satisfied; 

\\y^(j(t))l<  Myj(a))  l/fu  ;uid 

|>„0'wj|l  <  ( CO)  Vco  (^'-5) 

The  above  assumptions  imply  that  design  specifications  can  be  translated  into  the  bounds; 

o[Sg(jco)]  <  M^(co),  Vco  and  (C-6) 

ofSodoo)]  <  M^(0}),  Vco  (C-7) 

Since  there  is  a  conflict  between  achieving  these  two  goals  at  the  same  frequency; 

SJs)-hTJs)^  I  (C-8) 


there  is  a  trade-off  between  the  feedback  properties  governed  by  and  those  governed  by 
The  lack  of  directionality  properties  and  structure  in  the  systems  implies  that  this  trade-off  is 
essentially  the  same  as  that  in  the  SISO  case.  In  special  cases,  the  heuristic  rules  relating  the  open 
loop  gain  to  feedback  properties  can  be  extended  to  MIMO  systems.  Let  L(jco)  e  be  a 
transfer  function  evaluated  at  a  given  frequency  (o.  The  level  of  gain  experienced  by  an  input  to  the 
system  described  by  this  transfer  function  will  generally  depend  upon  the  direction  in  which  the 
input  lies.  Using  the  singular  value  decomposition; 

L=  (C-9) 

1^1 

If  the  open  loop  transfer  matrix  has  large  gains  in  all  directions  or  small  gains  in  all 
directions,  the  following  analogs  to  (C-1)  and  {C-2)  are  useful  in  assessing  the  feedback  properties 
of  the  system; 
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o[So(ja))] « 1 
^LJjco)]  »  1  <=>  and 

To(j0))==I  (C-10) 

and 

o[To(jO))] «  1 
«  1  <=>  and 

So(j(0)  =  I  (C-11) 

Multivariable  Design  SpeciHcations 

Suppose  that  sensors  in  some  loops  of  the  system  have  become  noisy  at  a  lower  frequency 
than  the  sensors  in  the  other  loops.  This  situation  can  be  mathematically  modeled  by  assuming  that 
the  sensor  noise  can  be  written  as  the  sum  of  a  stnictured  and  an  unstructured  component: 

n(s)  =  n^(s)  +  n^^(s)  (C-12) 

The  structured  noise  is  assumed  to  lie  in  a  k-dimensional  subspace  C  C”,  and  satisfies  the 
following  bound: 


\\n,(jo}]\\<MnJ(0).  Vcu 


(C-13) 


The  remaining  unstructured  component  of  the  noise  satisfies  the  bound: 

\\nJjcofi<M„ja))  Vco 

Suppose  that  the  ratios  M^Jco)  /  M^( o))  and  M„^(o})l  M^( (o)  are  as  depicted  in  Figure  C-1.  Let 
the  orthogonal  projections  onto  the  subspace  which  contain  the  structured  noise,  and  onto 

its  onhogonal  complement,  be  denoted  by  and  P/,  respectively.  Then  the  design 


specifications  for  this  system  should  satisfy  the  following  conditions: 

a[SQ(ja))]«\  (y  <  &»j 

(C-1 5) 

'^Jo(J(o)F,] «  1 

(1)  >  (02 

(C-1 6) 

^SoO<yiPv^] « i 

CO  ^  (U  3 

(C-17) 

o[To(  jco)] «  1 

<y  >  0)4 

(C-1 8) 

o[So(jo})]<  M^(a)) 

Vftj 

(C-19) 

o{T^(J(v)]  <Mj(CO) 

Vw 

(C-20) 
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The  above  specifications  could  be  satisfied  by  requiring  the  loop  gain  to  be  large  in  the  directions 
for  which  the  sensitivity  is  requind  to  be  small,  and  small  in  the  directions  for  which  the 
complementary  sensitivity  is  required  to  be  small.  To  describe  the  gain  in  different  directions,  the 
singular  value  decompxisition  of  the  loop  transfer  matrixL  is  partitioned  as: 

L  =  VyZyt/,"  +  V2l2^2  (C-2 1 ) 

wheie  Vy,{/y  ^  X;  =  d/flg[CT; . ,<7^]  and 

The  decomposition  requires  that  inputs  to  the  open  loop  system  lying  in  the  column  space  of  (/, 
are  amplified  by  the  level  of  gain  represented  by  the  singular  values  and  appear  as  outputs  in  the 

column  space  of  V,  This  partition  suggests  that  L  be  referred  to  as  consisting  of  a  higher  gain 
subsystem  VjZjU^  and  a  lower  gain  subsystem  V2Z2U2  ■  The  matrix  provides  a 

measure  of  the  alignment  of  the  column  spaces  of  Ui  and  V'',,  and  can  be  interpreted  as  a  measure 
of  interactions  between  the  high  and  low  gain  subsystems. 

Motivated  by  the  above  discussion,  con.sider  the  limiting  case  in  which  ojT;]  and  0(2,2] 
the  subsystem  decomposition  are  allowed  to  become  sufficiently  large  and  sufficiently  small, 
respectively.  The  following  theorem  shows  that  this  strategy  would  achieve  sensor  noise  rejection 
in  an  (n-/:)-dimensional  subspace  and  disturbance  rejection  in  a  /^-dimensional  subspace 
corresponding  to  the  directions  of  small  and  large  loop  gain,  respectively. 
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Theorem  C-l:  Let  detlulWA  ^  0.  Then,  if 


»  a] 


and 


o[^]«  ^^2^2] 


(C-22) 


(C-23) 


the  sensitivity  and  complementary  sensitivity  matrices  satisfy  the  following  conditions: 

Sensitivity  matrix:  ~  ^app]  '  and  ~  ^app]  ^ 

where 

So^=V2(vl'Uiy'v{‘ 

=■  ^2iy2^hy\^2^'\)^l'  -^^2^2 

-■■V,(v/'U2lV2^U2y\f+V2Vf 


Complementary  sensitivity  matrix: 


Tapp  =  V'l(yi"v,)‘'t^l" 

=  V, )■'((;, "^2)^2"+^^/' 


(C-25) 


The  above  theorem  allows  feedback  properties  to  be  approximated  when  open  loop  gain  is 
both  large  and  small  in  different  directions  at  the  same  frequency.  However,  one  must  take  caie  to 
correctly  identify  the  input  and  output  directions  of  the  high  and  low  gain  subsystems  with  the 
subspace  of  C”  in  which  the  signals  are  to  be  rejected. 

Next  consider  two  additional  cases:  when  the  gain  is  assumed  to  be  either  (a)  large  in  some 
directions  or  (b)  small  in  some  directions,  and  in  each  case,  no  assumption  is  placed  upon  gain  in 
orthogonal  directions,  the  following  two  theorems  show  that  there  exists  an  analogue  to  scalar  gain 
crossover  frequency  range  in  which  knowledge  of  gain  alone  is  insufficient  to  allow  approximation 
of  feedback  properties;  hence  it  is  necessary  to  consider  some  sort  of  phase  information. 
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Theorem  C~2.  Let  0  and  de/[W^t/2  +  X2]  0  .  Then,  if 

V/v,r 

^X,]  »  max\a^(u(wy'^l  -(ut^V2)£2(V2U2  -  -2^)‘‘(W 

5|((/,»F,)''(t/,''V2)X2(W2''(/2  + 


(C-26) 


then  the  sensitivity  and  complementary  sensitivity  matrices  satisfy  the  following  conditions; 


^  ^d  ^  (C-27) 

where  =  U2  (V^U2  +  ^2" 


and 


T  =  I  -S 
‘app  '  ‘^’^pp 


=  V,{ul'Viy‘L\'  +t,'2(V2»(.'2)'‘x2(v'2''t^2  +-^2)'’ 


V2 


Theorem  C.3.  Let  der[v2'^fy2]  ^  0  and  (:ter[/ +  f/f'v,Xi]  0.  Then,  if 


(C-28) 


(C-29) 


the  sensitivity  and  complementary  sensitivity  matrices  satisfy  the  following  conditions: 

5(^  -  S^pp] «  1  and  Sfr  -  Tppp]  «  1 

where  Sppp  s  / -Tppp  =  v!l  (V^V2)~'v^  +  (! +ui'v,Z^f'ul' 

Tapp  =  r,  f'ul'  (C-30) 
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Remark:  Theorems  C-3  and  Theorem  C-2  show  that  the  conditions: 

^£2) «  and  2[£2) »  5[y2''V2] 

and  the  conditions: 

ql^i]  >>  and  a[Li] 

fail  to  hold  at  the  same  frequencies.  Since  the  knowledge  of  the  open  loop  singular  value  gains 
does  not  piovide  sufficient  information  to  allow  feedback  properties  to  be  approximated,  additional 
information  about  the  phase  is  necessary.  This  area  requires  further  research. 

In  conclusion,  the  third  approach  by  Freudenberg  and  Looze  (1986),  (1990),  and  (1988) 
have  introduced  a  general  concept  of  the  relationship  of  the  gain,  phase  and  direction  of  the  open 
loop  transfer  matrix  with  the  closed-loop  sensitivity  and  complementary  sensitivity  matrices.  The 
rules  developed  in  each  of  the  three  theorems  can  serve  as  a  guide  for  design  and  robustness 
testing.  However,  more  direct  rules  still  need  to  be  developed  for  robustness  assessment  of  MIMO 
system.s. 
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